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_* i ABSTRACT 

We discuss some recent developments in the evaluation of the QCD corrections 
to Ap, their interpretation, an estimate of the theoretical error, and its effect on 
^ ■ electroweak physics. 
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The p parameter is conventionally defined as the ratio of the effective neutral and charged- 
current couplings at q 2 = 0. Although in the Standard Model (SM) the associated ra- 
diative correction Ap = 1 — p~ l is process-dependent, [] its fermionic component, (Ap)f, 
is universal. We recall that its contribution, pf = [1 — (Ap)/] -1 , is frequently separated 
out, as an overall renormalization factor, when neutral currents amplitudes are expressed 
in terms of G^, even when one considers amplitudes with q 2 ~ M§. Furthermore, (Ap)/ 
is related to the leading asymptotic corrections, for large m t , of the basic radiative cor- 
rections Ar, Af, and Ap Jl|@. For these reasons, it is highly desirable to evaluate (Ap)/ 
as accurately as possible, and to obtain an estimate of the theoretical error due to the 
unknown higher order corrections. In particular, the study of the QCD corrections to 
(Ap)/ has recently been the subject of considerable attention. 

Neglecting higher order electroweak effects ~ (a j 'n s 2 ) 2 \m 2 j 'M^) 2 , but retaining QCD 
corrections, (Ap)/ can be written as 

( A '>' = W [1+i « CDl ' (1) 

where m t is the pole mass, the first factor is the one- loop result ||, and Sqcd represents 
the relevant QCD correction. The 0(a s ) contribution was obtained some time ago by 
Djouadi and Verzegnassi || 

Sqcd = -^(vr 2 + 3)o(mt) + ... = -2.860a(m t ) + ... (2) 

where a(p) = a s (p)/ii. Attempts to go beyond this result have been carried out via two 
different methods: the dispersive approach, pioneered by Kniehl, Kuhn, and Stuart ||, 
and the direct examination of the relevant Feynman amplitudes at q 2 = 0. The dispersive 
approach and its comparison with other calculations is reviewed in the contribution of B.A. 
Kniehl to these Proceedings and in Ref . || . The present contribution only discusses recent 
developments in the Feynman diagram approach, their interpretation, and an estimate of 
the theoretical error. 

Very recently, Avdeev, Fleischer, Mikhailov, and Tarasov [[j] carried out a complete 
three-loop calculation of 0(aa 2 ). Their result, obtained in the limit m& — > 0, can be 
expressed to good accuracy as 

Sqcd = c x a(m t ) + c 2 a 2 (m t ) + ..., (3) 

where C\ = —2.860... is the Djouadi- Verzegnassi result, and 

c 2 = -21.271 + 1.786AT/ = -10.55... (4) 



1 See, for example, the second paper in Ref.[l]. 



In Eq. (H) Nf = 6 is the total number of quarks contributing via the vacuum polarization 
loops (5 massless quarks and the top). The contribution of the top quark is very different 
from the that of the massless quarks and has been split in two parts: one, corresponding 
to a "massless top" (hence the factor Nf = 6) and the remainder, which is included in 
the first term. As pointed out in Ref.[[jJ, it is also interesting to note that about 40% of 
C2 arises from the "anomaly-type diagrams", where the t^^t currents are attached to 
triangle diagrams linked by two virtual quarks ||. 

In order to discuss the evaluation of a s (/j) in these expressions, we recall that this 
parameter conventionally evolves with five active flavors for mj < \x < rrit, and with six 
for ft > m t . At // — m t , a s (//) is continuous through 0{a 2 s ). There is a very small disconti- 
nuity, — (25/72)a 3 (m t )0, but that occurs beyond the order of current (Ap)/ calculations 
and, moreover, is negligibly small. For the purpose of our discussion, we can therefore 
treat a s (/x) as being continuous at // — m t . From these observations it follows that a s (m t ) 
in Eq.(^j) should be evaluated conventionally, evolving a s (/i) from a s (M z ) with a 3- loop 
/3 function and five active flavors. As (Ap)j and m t are physical observables, they are 
//-independent. Therefore the same is true of Sqcd- However, because Eq.(|3]) is truncated 
in 0(af), its evaluation depends somewhat on the chosen scale. The //-dependence of 
Eq.(||) can be studied by using the simple relations 



a s (m t 






1-M 5) 
1-M 



71 



ln(- 

7r m t 



it m t 



+ ... 

+ ..., 



(5) 



Mf) 



S n t) 



,("t) 



where /3 1 n/ = — 1(11 — ^), the ellipsis stand for higher order terms, and the superscript 



nf in a s f and j3 1 f represents the number of active flavors in the evaluation of a s (//). 
Eq. (5) follows from the fact that a^ (p) and af^ (//) satisfy the RG equations with nf = 5 
and rif = 6, respectively, and from the above mentioned continuity at jjl = m t . Some 
authors employ af\ix) for \x < m t as well as for \x > m t . That choice, although the- 
oretically acceptable, is inconvenient because the experimental value of a s [Mz) should 
be identified with af\Mz) rather than a^(Mz). Therefore, it is natural to use the 
expression involving a^(fx) in Eq.(^) when m b < fi < m t , and that in terms of af\fx) for 
ji > m t . For clarity, the following related observations are relevant: i) The suggestion has 
been made by some theorists that, in replacing a s (m t ) by the r.h.s. of Eq.(|5|) with rif = 5 
and fi < m t , one should also change Nf from 6 to 5 in the //-independent coefficient c%. 
We see, however, that such procedure is theoretically inconsistent. While the shift from 
rif = 6 to rtf = 5 in Eq.(|J) is re-absorbed in a redefinition of a s (/t) in a /i-independent 
manner, a change of Nf in Eq.(||) for /i < m t makes Sqcd //-dependent in 0(a 2 s ). This 
contradicts the basic requirement that Sqcd should be //-independent through the order 



of the calculation. Or, putting this in a slightly different language, m t is //-independent 
and, unlike m t (//), it cannot be adjusted to absorb a /i-dependent contribution from Sqcd- 
ii) Some theorists evaluate a s (m t ) by evolving a^(/j.) from a s (M z ). As a s (M z ) should 
be identified with af'(M z ), we see from Eq.(|J) that such approach introduces a small 
but unnecessary error. 

In this report we follow the above observations: in the range mj, < // < m t , a s (/f) 
is evaluated by evolving a s (/x) from a s (M z ) with a three- loop j3 function^ and five ac- 
tive flavors. For definiteness, we take a s (M z ) = 0.118 and M z = 91.19GeV, and adjust 
A^ accordingly. In order to discuss the //-dependence of the truncated series or to opti- 
mize the perturbative expansions, we employ Eq.(||), choosing the five-flavors expression 
in the domain rn\, < H < vdt- As an illustration, for m t = 200GeV, Eq.(|3|) leads to 
Sqcd = —0.0961 — 0.0119 = —0.1080, so that the 0(a 2 s ) term implies an enhancement 
0.0119/0.0961 = 12.4% of the leading QCD result. If a s {m t ) in Eq.(|) is expressed in 
terms of a s (fi) and the resulting series truncated in 0(a 2 ), for 0.1 < fi/m t < 1 we find 
a variation < 5.8 x 10~ 3 , a bound that amounts to 5.4% of the total QCD correction or 
49% of the 0(a 2 s ) term. It is difficult, however, to make a precise statement of the error 
based on these considerations because the //-dependence of the truncated series depends 
on the chosen interval. For example, in this case the variation is significantly smaller in 
the interval 0.15 < fi/m t < 1 and significantly larger in the range 0.075 < ///m t < 1. 
Furthermore, the mildness of the /i-dependence of the first N terms of a QCD expansion 
over an interval is consistent with but does not imply the smallness of the higher order 
terms. The point is that, for example, the //-dependence induced by a significantly large 
0(af +1 ) term is of 0(a^ +2 ) and may be cancelled by sizable contributions of that and 
higher orders. The application of the Brodsky-Lepage-Mackenzie (BLM)|]lT|, Principle of 
Minimal Sensitivity (PMS)(T^|, and Fastest Apparent Convergence (FAC)|I3| methods to 
optimize the scale in Eq.(|3]) is discussed later on. 

On the other hand, we note that the expansion in Eq.([3]) involves rather large and 
increasing coefficients, a feature that frequently indicates significant higher order terms. 
Furthermore, arguments advanced in Ref. |l | suggest that there are at least two scales 
in Sqcd'- one, of 0(m t ), associated with the intrinsic corrections to the electroweak am- 
plitude, and another one, much smaller, related to contributions to the pole mass mt 
involving small gluon momenta. It is therefore a good idea to find alternative expressions 
for Sqcd that separate the two scales, and at the same time involve terms of 0(a 2 ) with 
coefficients of 0(1) rather than O(10). The advantage of this strategy is explained later 
on, when we discuss the theoretical error. A simple way of implementing this idea has 
been outlined in Ref.[|U|: one expresses first (Ap)j in terms of rht(m t ), the running MS 
mass evaluated at the pole mass, and then relates rht(mt) to mt by optimizing the expan- 



2 See, for example, Eq.(A.2) of Ref. @. 



sion for mt/rht{mt), which is known through 0(a 2 s ). Recently these arguments have been 
considerably refined |16f using the new results of Ref. 0. 

Calling \x t the solution of rh t (fi) = p and using Eq.(19) of Ref. [0], we have 



(Ap)j 



3£*llHt N , eMS 

Sv^vr^ [i + °' 



QCDh 



(6) 



5%% D = -0.19325 a(/2 t ) + 0.07111a 2 (^). 



(7) 



We note that the convergence pattern of Eq.(^) is very nice, with very small leading and 
next-to-leading coefficients. For this reason we will assume that the terms of 0(a^) and 
higher are negligible in Eq.(|6|) and evaluate Eq.(^) with fi t — > m t as this introduces only 
a small change of O(af). 

We can now express (Ap)/ in terms of rht(mt), by using the NLO expansion [ |T6| 1 



^ , = 1 + ^V*) + [35.96 - 2A5n f ]a 3 (m t ) 
m t (m t ) 3 

where n/ = 5 is the number of light flavors. Defining Aqcd by 

3G At m 4 2 (m < ) 
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we have 
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Neglecting terms of 0(af), we can combine Eqs.(|6|-|T0D in a single expansion: 

Aqcd = -0.19325 a(m t ) + Ca 2 (m t ). 



(10) 
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In Ref . |15| Eq. (|Tl|) was proposed before the results of Ref . ]7| became known, and it was 
argued, on the basis of convergence assumptions, that \C\ < 6. From Eqs.[7|-| we see that 
C = 16/3 + 0.071 = 5.40, consistent with the arguments of Ref. |L5[ -P| On the other hand, 
with C = 5.40 the two terms in Eq.flTTD nearly cancel and the error estimate is unneces- 
sarily large. It is clearly better to insert in Eq. ([T0|) the separate expansions in Eqs.(0)and 
(§), as the magnitude of the last terms in these expressions is significantly smaller than 
in Eq. (|TTD . In the case of Eq.(|^) we retain the relatively large 0(a 3 ) contribution in order 



3 In Ref. jl5| C was also estimated to be w +3 by optimization arguments, but the more conservative 
value C = 0^6 was employed in the final analysis. 



to control the scale of the leading term. As its coefficient, 23.71, is rather large, we apply 



the BLM optimization procedure JTT| to Eq.flg), and obtain 

^ r = 1 + -a 2 (0.252m t ) - 4.47a 3 (0.252m t ) (BLM). (12) 

m t {m t ) 3 

As an illustration, for m t = 200GeV, Eqs.(||) and ( |i~2"|) give fj, t /rht(m t ) = 1.00391 and 
1.00423, respectively. The optimization of Eq.(|) using the PMSffT2| and the FAC p|| 
methods is very close numerically to Eq. ([12]) (the difference is < 5 x 10 -6 ). 

Values for Sqq d (Eq.(^) with \i t — > m t ) and Aqcd (evaluated via Eqs.(|T^), (0) and 
(0)) are shown in Table 1. We see that these are indeed small effects. In particular, 
Aqcd = (2 — 3) x 10~ 3 , depending on m t , a remarkably small correction. Comparing 
Eq.(l) and (@), we find 

1 + 5 QCD = {^) 2 '[I + Aqcd]. (13) 



For m t /rht(mt) we have the well-known expansion [J17 



m t -, , 4 ,.. 2 



1 + -a(m t ) + K a 2 (m t ), (14) 



m t (m t ) 3 

In the limit of neglecting the masses of the first five flavors, K is given by PIJ18 



tf = 16.0065-n/1.0414 + 0.1036, (15) 

where the first term corresponds to the quenched approximation, and the second and 
third are the vacuum polarization contributions of the masseless quarks (n/ = 5) and the 
top quark, respectively. As expected, the latter is much smaller. Refs.f7|,|]l7| include the 
contribution of all flavors to the vacuum polarization and, therefore, in those calculations 
rht([i) evolves with six active quarks in the (3 and 7 coefficients. The same is true of 
m t (/u t ) = fi t in Eq.(|) and rht(m t ) in Eq.([L4]). One can equivalently express these relations 
in terms of rh t (fit) and rht(mt) evaluated with five active flavors 0. Through terms of 
0(a 2 s ) this is simply done by "decoupling", i.e. subtracting the small top contribution, 
0.1036, from K and 2 x 0.1036 w 0.21 from the coefficient of a 2 in Eq.(0)fj These 
changes are very small and, moreover, they are not necessary for our purposes. As our 
aim is to combine Eqs.(@) and fll4D , both evaluated with six flavors, we have defined in 
the same way fi t and m t {m t ) in Eq.(|[) et seq., without decoupling the small top quark 
contribution. 



4 It is worth emphasizing that even this small decoupling subtraction should not be made in discussing 
the /i-dependence of Eq.(l) because, unlike TO t (/i), mt is \x- independent. 



When Eq.(|T3D is inserted in Eg. (pT3|) , it induces a contribution ~ —0.11 to Sqcd while, 
as we saw before, Aq C d ~ 0.002. Thus, we find the intriguing result that, when (Ap)f is 
expressed in terms of m t , the correction is almost entirely contained in (rh t (m t )/m t ) 2 , a 
pure QCD effect that can be studied in isolation from electroweak physics. Once this is 
recognized, it becomes clear that the magnitude and error of Sqcd are largely controlled 
by the value of (rn t (m t )/m t ) 2 and the accuracy within which it can be calculated. As K 
is quite large (10.90 for n/ = 5), it is natural to apply the three well-known approaches 
H|-||13|| to optimize the scale in Eg. (|14D . Given an 0(a 2 ) expansion of the form S = 



1 + Aa(m t ) + (B — Crif)a 2 (mt), where A, B, and C are numbers, the application of these 
methods lead to the optimized expansions 

33 

S = 1 + A a(p*) + [B- —C}a 2 (p*) (BLM), (16) 



S=l+Aa(/i tt )-^aV) (PMS), (17) 

2pi 



S = 1 + A a(p***) (FAC), (18) 

where (3 2 = -[51 - 19n//3]/4, /i* = e- 3C/A m t , /i*** = e {B ~ Cn f )/Af3l m t , and //** = 
e (/V 2 /9i)//**, In deriving these expressions we have assumed that a s evolves with rif 
active flavors between m t and the optimized scales. To obtain Eq.([T7|) we have employed 
the two- loop RG differential operator to find out, to good accuracy, the stationary point. 
This can be checked by inserting Eq.([|) in the r.h.s. of Eq.flHD and evaluating numerically 
the series, truncated in 0(a 2 ), as a function of \i. Applying Eqs. fli6| - |l8"| ) to Eq. (|14|) , 



,; '' ! ' A -a(fi*) - 1.07aV) (BLM), (19) 



m t (m t ) 3 



m < 1 i ^ _ , :,• , ,< .-iA ,2/,,** 



m t (m t ) 3 



l + -a(/i")-0.84a 2 (/i**) (PMS), (20) 



,/;; l + ^a(fi***) (FAC), (21) 



rh t (m t ) 3 

where /i* = 0.0960m i5 /j,** = 0.1005m t , fi*** = 0.1185mt- In these expressions a s is 
evaluated with n/ = 5. As the optimization scales are in the m& — m t range, this 
is consistent with the discussion after Eq.(^). For m t = 200GeV, Eqs.( |T9"[ - PTD give 
">TT't/ r > >> 't( m 't) — 1-06304, 1.06303, and 1.06297, respectively. For m t = 174GeV, the cor- 



responding values are 1.06477, 1.06477, 1.06470. Finally, for m t = 130GeV, we have 
1.06875, 1.06875, 1.06867. Thus, the three approaches give expansions with similar scales, 
coefficients of 0(1), and remarkably close values. In contrast, the expansion in Eq. (|14|) , 
which involves a large second order coefficient of O(10), gives 1.0571, 1.0584, 1.0614 for 
m t = (200, 174, 130)GeV, which are (0.6 — 0.7)% smaller. One must conclude that the 
coefficients of the unknown terms of 0(a 3 ) and higher in Eq.(|14]) and or Eqs. ( p!9] - pl| ) are 
large. For instance, if Eq.(|l9|) were exact, the coefficient of the a 3 (m t ) and a 4 (mt) terms 
in Eq.(|l~4"D would be m 104 and 1,041, respectively.^ In the following we employ the 
optimized expression for mt/rht(mt), which, for definiteness, we identify with the BLM 
expansion (Eq.([T9|)). The advantage of this procedure will be explained later on, when 
we discuss the theoretical error. 

Table 2 displays the values of Sqcd obtained from Eq. flT3"|) , using Eq.(|T9"D and our 
previous determination of Aqcd (Eq. ([T0|) ), and compares them with those derived from 
Eqs.(0, |J). In order to show the effect of the higher order contributions (H.O.C), we 
also exhibit the fractional enhancement of the total QCD correction over the conven- 
tional O(a) result (Eq.(|)). From Table 2 we see that, for 130GeV< m t <220GeV, the 
evaluation of |<5<2od| from Eq.(|I~3"D leads to a (18 — 20)% enhancement over Eq.(fJ) and is 
(5.2 — 6.6) x 10~ 3 larger than the results from Eq.(^). As a percentage of the total QCD 
correction this difference amounts to ~ 5%, a reasonably small effect. However, the last 
two columns in Table 2 show that the H.O.C. are ~ 45% larger in the evaluation based on 
Eq . (|i~3"D . Most of this divergence can be traced to different ways of treating the dominant 
factor {rht(m t ) / m t ) 2 in Eq. (|T3"D . For example, it is amusing to note that the expansion 
{m t {m t ) I m t ) 2 — 1 = — (8/3)a(m 4 ) — {2K — 16/3)a 2 (mt), obtained from Eq . (|i~4"D , is quite 
close numerically to Eqs.(|, §. For m t = (200, 174, 130)GeV it gives -0.1081,-0.1106, 
-0.1160, to be compared with -0.1080, -0.1102, 0.1154 from Eqs.(|, § (see Table 2). This 
conforms with our observation that the bulk of the QCD correction to (Ap)j is contained 
in (rht(m t )/m t ) 2 , but at the same time shows the main source of the difference between 
the two calculations. In fact, in the evaluation of Eq . (|l"3"D we have employed Eq.(|i~9"D, the 
optimized version of Eq . (|i~4"D , rather than the expansion mentioned above. 

As C2 in Eq.(||) is O(10), in analogy with our discussion of Eq. (]T4]) , we may directly 
attempt to optimize Eqs.(^|^). Applying Eqs. (|T9| - pTD to Eqs.(|3]}4|) with Nf — nf + 1 = 6, 
we have 

S QCD = -2.860 a(0.154m t ) + 9.99 a(0.154m t ) (Eq.(3); BML), (22) 

S QCD = -2.860 a(0.324m t ) + 1.80 a 2 (0.324m t ) (Eq.(3); PMS), (23) 



5 For recent applications of the BLM, PMS, and FAC approaches to estimate higher order coefficients 
in other cases, see Ref . |Ol , [Eof . 



S QCD = -2.860 a(0.382m 4 ) (Eq.(3); FAC). (24) 



Once more, in Eqs. (|22| - |24|) a s is evaluated with 5 active flavors, in accordance with the 
discussion after Eq.(@).Q For m t = 200GeV, Eqs.@^2|) give 0.1084, 0.11045, 0.11038, 
respectively. The difference between Eqs. (23|,24]) and Eq. ([T3|) is now only « 2.9 x 10 -3 , 
which is ~ 45% smaller than that between Eq.(|]) and Eq. (|T3"D . On the other hand, 
although the a 2 coefficient in Eq. (p3|) is roughly of 0(1), that in Eq.(|22|) is O(10). It has 
been pointed out that in many NLO QCD calculations, it is a good approximation to 
retain only the leading term, evaluated at the BLM scale ||L1||,||14||. From Eq.(p2[) we see 



that this is not the case in the expansion of Eqs.(||[|). Thus, in contrast to Eqs. ([l9| - pTT ), 
when applied to Eqs.(|3],fD the three optimization procedures do not uniformly lead to 
coefficients of 0(1) and similar scales. 

A frequently used method to estimate the theoretical error of QCD expansions is to 
consider the magnitude of the last included or known term Plfl . By this criterion the 
dominant error in Eq.flTB]) is contained in the first factor. We employ Eq.(|T9|) as it is the 
optimized expansion with largest second order term. The corresponding error is 

S f^V « ±] » im f { C - ±1.77aV). (25) 

\ m t ) (m t /m t (m t )y 

For m t = (200, 174, 130)GeV, this amounts to ±(4.3,4.5,5.1) x 10~ 3 . The calculation 
of the small correction Aqqd involves in turn two uncertainties. One, associated with 
(Ht/rhtirrit)) 2 , can be estimated from Eq.([T2|) as 

S ((/i t /m t (m t )) 2 ) w ±8.9a 3 (0.252m t ), (26) 

and leads to ±(6, 7, 8) x 10~ 4 . The other, involving Sq^ d , is given by 

<K<& « ±0.071a 2 (m t ) (27) 



and amounts to ±(8,8,9) x 10~ 5 . In Eq.(|l3]) these two errors are decreased by a factor 
(rht(mt)/mt) 2 ~ 0.88. Adding the three uncertainties linearly we get an overall error 
estimate for Eq.(^) due to higher order corrections: 

5(6 QC d) = ±(4.9, 5.2, 5.9) x lO" 3 (Eq. (13)) , (28) 

for m t = 200, 174, 130GeV, respectively. It could be argued that the smallness of the 
leading and NLO coefficients in Eq.(^) is fortuitous. However, in order to lead to a 



6 In Ref.|7J the FAC scale is reported as being 0.348mt. The apparent difference with Eq.(p4[) is due 
to the fact that a B is employed in that paper. In fact, both results are equivalent through 0(a 2 s ). 



contribution m 2.5 x 10 -3 for m t = 200GeV, which would modify our error estimate by 
50%, the coefficient of a 3 (m t ) should be ~ 66, or 930 times the coefficient of a 2 (m t ). 
If the same "last term" criterion were applied to Eq.(^), the error estimate would be 
S(Sqcd) ~ ±C2a 2 (m 4 ), which amounts to ±(1.2, 1.2, 1.3) x 10~ 2 . This is 2.2-2.4 times as 
large as in Eq . (|13"D . The consideration of the optimized expansions is very ambiguous in 
the case of Eqs. ([3], |j); while the PMS method (Eq.(pB"D) leads to a small error estimate, 
the BLM approach (Eq.(^)) gives a very large uncertainty p» ±(2.0 — 2.4) x 10~ 2 . This 
ambiguity may perhaps be related to the observation that the expansion in Eq.([3]) involves 
more than one scale. 

We note that the two different evaluations of Sqcd, by means of Eqs. (|||) and Eqs. ( |i~3]JiO| JK^ 



respectively, coincide through 0(a 2 {m t )). The numerical difference between the two 
means that at least in one of these two calculations there are significant contributions 
of 0(a 3 ) and higher. Although the present author prefers the latter approach, on the 
grounds that it involves expansions with coefficients of 0(1) and leads by the "last term 
criterion" to a significantly smaller error estimate, it seems impossible to rigorously decide 
at present which is more accurate. This could be clarified, in principle, by evaluating the 
0(a 3 ) terms in Eq.(|3]) and or Eq.([b|). Prospects for achieving this, however, appear to 
be quite remote [Tj|. On the other hand, we believe that Eq. (|28D is a reasonable estimate 
of the theoretical error due to unknown higher order corrections. In particular, we note 
that Eq. (|28|) is also very close to the difference between the two Sqcd evaluations, which 



can also be used as an estimate of the theoretical error [jTg, and to the scale variation 
of Eq.(|3]) in the interval 0.1 < fi/m t < 1. Finally, one should remember that there is 
at present an additional 5% error in Sqcd associated with the ±0.006 uncertainty in 
a s (M z ). This leads to an additional contribution of ±(5.8, 5.9, 6.2) x 10~ 3 to S(Sqcd), for 
m t = (200, 174, 130)GeV. In contrast to the theoretical error in Eq.(^), it is likely that 
this uncertainty can be decreased in the near future. 

As pointed out in Ref.|16[, the values for Sqcd given in Table 2 can accurately be 
represented by simple empirical formulae. We find that the calculation based on Eq.(13) 
(third column of Table 2) and the error estimate of Eq. (28) due to higher order corrections 
can be conveniently expressed as 

Sqcd = -2.860a(fm t ), (29) 



£ = 0.32118$° (Eqs.(13,28)), (30) 



while Eq.(3) (second column of Table 2) corresponds [0| to £ = 0.444. We emphasize 



that Eqs. (29,30) are not the result of a FAC optimization. They are simply heuristic 
formulae that reproduce the values in the Tables with errors of at most ±1 x 10 -4 for 
£ = 0.321 and £ = 0.431, and at most 3 x 10~ 4 for £ = 0.248. We also emphasize that in 

9 



these expressions a s is evaluated with five active flavors in the manner explained before. 
We see that Eq.(30) is somewhat smaller than the effective scale associated with Eq.(3); 
however, both evaluations are roughly consistent within our theoretical error estimate. 
By a numerical coincidence, the central value in Eq.(30) is very close to the PMS scale in 
Eq.(23); however, Eqs.(29) and (23) differ somewhat because of the presence of the 0(a 2 ) 
term in the latter. 

In summary, in this report we have emphasized a number of results and observations: 
i) Working in the NLO approximation, we have pointed out that Aq C d, the QCD correc- 
tion when (Ap)f is expressed in terms of rh t (m t ), is remarkably small, ~ (2 — 3) x 10 -3 . 
This means that precision electroweak physics essentially measures this parameter, an 
important input for GUTs studies. The evaluation of Aqcd and <w«d shows that rht(m t ) 
and fit = Tnt(fJ>t) are very good expansion parameters as they absorb the bulk of the 
QCD corrections |p2[. Arguments for these desirable properties from effective field the- 



ory have been given by Peris ]23|]. If the smallness of Aqcd persists beyond the NLO 
approximation, it also means that Sqcd, the QCD correction when (Ap)f is expressed in 
terms of mt, is almost entirely contained in (rht(mt) / 'mt) 2 , a pure QCD effect that can be 
studied in isolation from electroweak physics, ii) We have also pointed out that when the 
BLM, PMS and FAC optimization methods are applied to the expansion for m t /rht(mt), 
they lead to similar scales, coefficients of 0(1) and remarkably close values. Combining 
these optimized expansions with Aqcd one obtains an evaluation of Sqcd which shows 
an (18 — 20)% enhancement over the two-loop calculation, depending on mt- Using the 
magnitude of the last included terms as an estimate of the theoretical error due to higher 
order corrections, this analysis leads to S(Sqcd) ~ ±(5 — 6) x 1CT 3 , depending on m t . This 
estimate seems reasonable, as it is also close to the difference between the Sqcd evaluation 
mentioned above and that obtained from the expansion proposed by Avdeev et al.. If the 
theoretical error is combined quadratically with that arising from Sa s one obtaines an 
overall uncertainty in 5 QC d of ±(7.6,7.9,8.6) x 1(T 3 for m t = (200, 174, 130)GeV. For 
m t = 200GeV and M H = 300GeV, this induces errors ±9.5 x 10~ 5 in Ap, ±3.1 x 10~ 4 
in Ar, ±3.2 x 10~ 5 in the MS parameter sin 2 6 W (M Z ), and ±840MeV and ±5.5MeV in 
the predicted masses of m t and M w , respectively. We also recall that these results scale 
approximately as m 2 , so that they are actually smaller for m t < 200GeV. Thus the effect 
of these errors in electroweak physics is rather mild. It should also be remembered that 
the concept of pole mass has an intrinsic uncertainty ~ Aqcd and, for the top quark, this 



may amount to 200-300MeV [24 
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mt(GeV) 


io 3 fe 


10 3 A QCD 


130 


-6.80 


2.92 


150 


-6.67 


2.59 


174 


-6.53 


2.28 


200 


-6.41 


2.01 


220 


-6.33 


1.84 



Table 1. The corrections 6q^q d and Aqcd- The first one is given by Eq.(7) with \Xt — > rrit, 
while the second is obtained from Eq.(10), with fi t /m t (m t ) evaluated according to Eq.(12) 
(a s (M z ) = 0.118 is employed). 



m t 
(GeV) 


Sqcd 
(Eq.(3)) 


8qcd 
(Eq.(13)) 


H.O.C. 

(Eq.(3))% 


H.O.C. 

(Eq.(13))% 


130 


-0.1154 


-0.1220 


13.2 


19.6 


150 


-0.1128 


-0.1189 


12.9 


19.0 


174 


-0.1102 


-0.1160 


12.6 


18.5 


200 


-0.1080 


-0.1133 


12.4 


18.0 


220 


-0.1064 


-0.1116 


12.2 


17.7 



Table 2. Comparison of two determinations of Sqcd- The second column is based on 
Eq.(3) [6]. The third column is based on Eq.(13) with m t /rh t (m t ) obtained from Eq.(19) 
and Aqcd evaluated according to Table 1. The fourth and fifth columns give the fractional 
enhancement over the conventional 0(a s ) result (Eq.(2)) due to the inclusion of higher 
order contributions (H.O.C). 
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